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ANALYSIS OF THE INJECTION OF A HEATED TURBULENT JET 

INTO A CROSS FLOW* 

By James F. Campbell and Joseph A. Schetz^ 

Langley Research Center 

SUMMARY 

An investigation has been undertaken to develop a theoretical model of the incom- 
pressible jet injection process. The discharge of a turbulent jet into a cross flow was 
mathematically modeled by using an integral method which accounts for natural fluid 
mechanisms such as turbulence, entrainment, buoyancy, and heat transfer. The analyt- 
ical results are supported by experimental data and demonstrate the usefulness of the 
theory for estimating the trajectory and flow properties of the jet for a variety of injec- 
tion conditions. The capability of predicting jet flow properties, as well as two- and 
three-dimensional jet paths, was enhanced by obtaining the jet cross-sectional area dur- 
ing the solution of the conservation equations (a number of previous studies assume a 
specific growth for the area). Realistic estimates of temperature in the jet fluid were 
acquired by accounting for heat losses in the jet flow due to forced convection and to 
entrainment of free- stream fluid into the jet. 

INTRODUCTION 

One of the more basic processes in fluid mechanics is the mixing that occurs during 
the interaction of two intersecting streams of fluid. A large portion of such interaction 
processes exists where one stream of fluid has a much smaller mass (or volume) flow 
than the other stream, such as the case of a plume issuing from a smoke stack or a fluid 
injecting into a boundary layer on a vehicle traversing the atmosphere. This class of 
interaction problems, categorized as fluid or jet injection processes, is usually asso- 
ciated with a nearby solid boundary or surface. These injection processes are compli- 
cated by the fact that in their natural state they are almost invariably turbulent. 


*The information presented herein was included in a thesis submitted by James F . 
Campbell entitled "Analysis of the Injection of a Heated, Turbulent Jet into a Moving 
Mainstream, With Emphasis on a Thermal Discharge in a Waterway," in partial fulfill- 
ment of the requirements for the degree of Doctor of Philosophy in Aerospace Engineering, 
Virginia Polytechnic Institute and State University, Blacksburg, Virginia, December 1972. 

^Virginia Polytechnic Institute and State University. 
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Details of the jet injection process have been obtained by a number of experimenters 
(refs. 1 to 13) and reveal the very complex velocity and temperature fields that exist due 
to the interaction of the jet flow with the cross flow. Theoretical attempts (refs. 14 to 
22) to model this interaction process have usually considered the injected fluid as a tur- 
bulent momentum jet having a two-dimensional path. Since a large number of physical 
problems are concerned with jets having three-dimensional paths, it is desirable to have 
a theory that is capable of calculating three-dimensional jet trajectories as well as one 
which allows estimates of jet flow properties to be made. 

Accordingly, the present study was initiated to develop a theoretical method which 
accounts for the pertinent fluid-mechanic and heat-transfer aspects of heated turbulent 
jets discharging into a cross flow. The injection process is assumed to be incompressi- 
ble where the jet flow injects into a cross flow of semi-infinite extent and to have a spa- 
tially dependent velocity field. Emphasis has been given to determining the effects of the 
ambient flow variables and injection conditions on the jet path and to examining how the 
flow properties of and heat loss from the jet vary along the path. 

The general development of the theoretical model is presented in the main body of 
the report along with comparisons with other analytical methods and with experimental 
data acquired from a number of investigations. Details of the theory are contained in the 
appendixes . 


SYMBOLS 

A cross-sectional area of jet control volume 

Aj^ area of jet orifice, 

C effective jet circumference, sJAttA 

Cj) drag coefficient 

Cp pressure coefficient 

Cp specific heat at constant pressure of jet fluid 

D drag on jet flow due to blockage of free- stream flow 

d effective jet diameter, sJaA/t! 

dj diameter of jet orifice 
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differential length of jet control volume 


elemental area 
elemental volume 

entrained mass flow per unit length of jet 
entrainment coefficient 
unit vector normal to elemental area d6 
unit vectors in natural coordinate system 
unit vectors in Cartesian coordinate system 
pressure force 

body force acting on elemental volume do 
centrifugal force per unit volume 
buoyancy force per unit volume 
gravity acceleration constant 
average film heat -transfer coefficient 
width of jet control volume 

nondimensional entrainment parameter, see equation (B2) 

finite -difference grid point 

thermal conductivity of jet fluid 

mass of jet fluid in control volume 

free- stream mass entrained into control volume 



%r,i 


N. 


Nu,d 


N 


Pr 


N 


Re,d 


N 

P 

P 

Poc 

Q 

q 

R 


Re,di 


Froude number of jet flow at injection point, 
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Nusselt number based on effective jet diameter, Hd/k 
Prandtl number 

Reynolds number based on effective jet diameter, see equation (47) 

Reynolds number of jet flow at injection point, V^d^y/j^. 

local static pressure around perimeter of jet cross section 

average static pressure in jet flow 

free- stream static pressure 

rate of heat flow from jet control volume 

average dynamic pressure of jet flow, 

average dynamic pressure of free-stream flow, p V 2/2 

^ ^ 00 / " 

radius of curvature of jet trajectory 


r position vector from injection point to a point on the jet trajectory, 

see figure 32 

®c effective cylindrical area of jet control volume, TrdAs 

Sref reference area of jet control volume 

s,n,t natural coordinate system attached to jet trajectory, see figure 32 

T average temperature of jet fluid, see equation (4) 

Too average temperature of free-stream fluid 

t time, see equations (1) and (9) 
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u,w 

V 

V 


Ve 

Veff 

V,o 

x,y,z 

yo 


a,^,y 


As 

AT 

V 
9 
A 

V 


P 

P 


direction cosines, cos a and cos (3, respectively 
average velocity of jet flow, see equation (4) 


velocity vector 

entrainment velocity, see equation (B5) 


effective velocity ratio, 



average free-stream velocity 


Cartesian (inertial) coordinate system, see figure 2 

point above injection surface where laminar jet flow is assumed to begin its 
turbulent growth 

inclination of jet axis with respect to X-, Y-, and Z-axis, respectively, 
see figure 32 

infinitesimal length of jet control volume 
temperature difference, T - 
defines expression in equation (A22) 
angular orientation, see figure 8 
defines expression in equation (A23) 
coefficient of viscosity 
kinematic viscosity 

average density of jet fluid, see equation (4) 
average density of free-stream fluid 
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surface stress tensor 


T 

Tq torsion of jet trajectory, see equation (A9) 

constants in expression for see equation (46) 

Subscripts: 

i conditions at point of jet injection 

I local value of jet flow property 

nom nominal 

s,n,t conditions in s-, n-, and t-direction, respectively 

jet flow conditions in control volume before and after entrainment 
Cl)>(D>(D surfaces of jet control volume, see figure 2 
°° conditions associated with the free-stream flow 

Arrow over symbol indicates vector. 

REVIEW OF LITERATURE 

There are a variety of problem areas in incompressible flow situations that would 
benefit from understanding the basic fluid injection process. In aeronautical disciplines 
some of these areas are boundary-layer control, jet flap technology, fuel injection, aero- 
acoustics, and V/STOL aerodynamics. In the environmental sciences there is the dis- 
charge of pollutants into the atmosphere or into various water systems. Margason and 
Fearn's reference list (ref. 1) on previous experimental investigations of incompressible 
jets injecting into a cross flow covers some of these research areas and, thus, provides 
3.11 3.ppropri3t6 starting point for tho prosont review of relevant literature, 

The complex interaction that takes place after a fluid is injected into a moving 
stream results in a three-dimensional flow field, even though the jet may be following a 
two-dimensional path. (See fig. 1.) This, of course, complicates the job of experimen- 
tally measuring the complete flow field. Gordier (ref. 2), Margason (ref. 3), and Flatten 
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and Keffer (ref. 4) concentrated their efforts on meas- 
uring the trajectory of^the jet for a wide range of injec- 
tion angles and velocity ratios Vi/V^o, the trajec- 
tory being one of the easiest jet properties to measure. 

(See fig. 2.) Jordinson (ref. 5), Keffer and Baines 
(ref. 6), Ramsey (ref. 7), Kamotani and Greber (ref. 8), 
and Fricke, Wooler, and Ziegler (ref. 9) measured 
details of the interaction process and revealed the 
very complex flow field that exists due to the pressures 
and shear stresses in and around the jet flow. Very 
good descriptions of how the jet flow distorts and 

develops under the influence of the free-stream flow 

, Figure 1.- Diagram or tne inrerac- 

and body forces are given by Keffer and Baines (rel. bj, resulting from jet injec- 

Ramsey (ref. 7), Abramovich (ref. 10), and Keffer tion into a cross flow (afM;er 

Abramovich (ref. 10)). 

(ref. 11). 

The most detailed measurements of the flow field appear to be those of Keffer and 
Baines (ref. 6) and Kamotani and Greber (ref. 8). In particular, Kamotani and Greber 
examine the structure of the rotational velocity field which results because of the shear- 
ing action between the free-stream flow and the edge of the jet flow. The rotational veloc- 
ity field, usually interpreted to be a pair of counterrotating vortices, is an interesting as 
well as important facet of the injection process; it affects the path of the jet as well as the 
mechanisms that govern entrainment. In both reference 6 and 8, the data were reduced 
so that a measure of mass flux in the jet could be obtained, which in turn resulted in esti- 
mates of the rate that free-stream fluid is entrained into the jet structure. In addition to 
investigating the characteristics of the velocity field, Ramsey (ref. 7) and Kamotani and 
Greber (ref. 8) measured the temperature field resulting from the injection of a heated 

jet. 

It is important to note the limiting conditions of a jet injected at an arbitrary angle 
into a cross flow. If the angle between the jet axis and the direction of the free-stream 
velocity a goes to zero for a given value of Vj^/V^, the condition exists of a jet in a 
coflowing stream, that is, where the jet flow is parallel to the free stream. If the jet 
velocity becomes very large (Vj/V^ - «°) for a given a, the situation approaches that 
of a free jet. An indication of the magnitude of entrainment for a jet in a coflowing 
stream and for a free jet has been provided by Morton (ref. 12) and Ricou and Spalding 
(ref. 13), respectively. 

In reviewing some of the theoretical methods available for modeling an injection 
process, it should be emphasized that a method is desired here which allows the injected 
fluid to be followed from its point of discharge to some point downstream where complete 


Rotational 
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mixing takes place. At a minimum, from the method one should be able to estimate the 
three-dimensional (3-D) path of the injectant, to show how the jet size varies as it pro- 
ceeds downstream, and to include basic heat-transfer mechanisms which allow the jet 
temperature to be predicted. The most advanced method would result in a complete and 
detailed description of the flow field resulting from an injection into a moving stream. 
This description could be accomplished by solving the full three-dimensional, turbulent, 
Navier-Stokes equations in an Eulerian framework, which requires the specification of 
the eddy viscosity field. (See ref. 14, for example.) 

In order to avoid the complexities inherent with this approach, many studies have 
tried to theoretically model the gross features of the injection process by describing the 
fluid motion of the jet from the point of discharge in a Lagrangian framework. This pro- 
cedure allows an estimate of jet properties to be obtained if the appropriate forces acting 
on the jet flow are accounted for. Since the jet path is the most obvious of the jet proper- 
ties, it is natural that early attempts were concerned only with obtaining estimates of the 
trajectory. Abramovich (ref. 10), for example, obtained the trajectory of a jet which had 
a circular cross section at the injection point by balancing the centrifugal and blockage 
forces perpendicular to the trajectory. His basic argument was that the blockage effect 
of the jet flow on the free-stream flow could be approximated by assuming that the jet 
flow acts as a ''solid” body inclined at some angle to the free stream. He accounted for 
the deformation of the jet cross section by assuming the shape to be elliptical and by 
specifying a growth rate for the cross-sectional area, which is necessary if only one 
force equation is used to obtain a solution for the trajectory. One of the serious draw- 
backs of this method is the assumption that the component of jet momentum perpendicular 
to the direction of the free-stream flow remains constant. This assumption was relaxed 
by Schetz and Billig (ref. 15). 

Two other forces acting perpendicular to the jet trajectory and which help govern 
the development of the jet flow are (1) buoyancy force resulting from a difference in den- 
sity between the jet and free-stream fluids and (2) entrainment force resulting because of 
the free-stream fluid that is drawn into the jet structure. Theoretical trajectories were 
obtained by Reilly (ref. 16) and Campbell and Schetz (ref. 17) using procedures similar to 
those previously described but also accounting for the entrainment phenomenon; Campbell 
and Schetz, in addition, included the buoyancy force in their model. 

Since all these previous works utilized an assumed area growth based on experi- 
mental data obtained in the proximity of the injection point, they are not suitable for pro- 
viding realistic trajectory information farther downstream. This deficiency can be 
avoided if, instead of assuming an area growth rate, a momentum conservation equation 
in the direction of the jet path is used; this equation is used in addition to the conserva- 
tion equation in the direction perpendicular to the trajectory. Wooler, Burghart, and 
Gallagher (ref. 18), Hoult, Fay, and Forney (ref. 19), and Hoult and Weil (ref. 20) used 
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this procedure, solving force ecfuations normal and parallel to the jet path simultaneously 
to obtain a solution for the trajectory. An added advantage of using these two force equa 
tions is that, if all the appropriate forces are accounted for, the solution procedure allows 
the jet flow properties to be estimated, as was done by Hirst (ref. 21) and Campbell and 
Schetz (ref. 22). The forces accounted for in a number of theoretical investigations are 

summarized in table I. 

TABLE I.- CAPABILITIES OF THEORIES USING INTEGRAL TECHNIQUES 


■ 

Investigation 

Arbitrary 

injection 

velocity 

Arbitrary 

injection 

angle 

Forces parallel and perpendicular 
to jet trajectory 

Computed 
jet flow 
properties 

3-D 

trajectory 

Blockage 

Entrainment 

Buoyancy 

P, A,V 

T 

Abramovich (ref. 10) 

y 

y 

y 






Schetz fii Blllig (ref. 15) 

s/ 

y 

y 






Reilly (ref. 16) 

y 


y 

y 





Campbell & Schetz (ref. 17) 

v/ 


y 

y 

y 




Wooler, Burghart, & Gallagher 

y 


y 

y 





(ref. 18) 









Hoult, Fay, & Forney (ref. 19) 

y 



y 

y 




Hoult & Weil (ref. 20) 

y 



y 

y 




Hirst (ref - 21) 

y 

y 

y 

y 

y 

y 


^y 

Campbell & Schetz (ref. 22) 

y 

y 

y 

y 

y 

y 

y 


Present 

y 

y 

y 

y 

y 

y 

y 

y 


®^No 3'D trajectory results given. 


The theories discussed thus far have been concerned with predicting the flow char- 
acteristics of a jet following a two-dimensional (2-D) path. This type of path occurs when 
the radius -of -curvature vectors associated with the trajectory lie in one plane, the tra- 
jectory, of course, being a curve in that plane. If the injection and free-stream velocity 
vectors shown in figure 2 are considered to form an "injection plane," then a two- 
dimensional trajectory will always result when this plane is oriented vertically, that is, 
alined with the gravity vector. When the injection plane is rolled away from the vertical, 
a three-dimensional trajectory may occur depending on the buoyancy of the jet flow. If 
the buoyancy force is absent, or is small with respect to the jet momentum, then a two- 
dimensional path will result. A larger buoyancy force, however, will cause the jet to 
bend out of the injection plane, and the result will be a three-dimensional trajectory. 

The radius-of- curvature vectors associated with this type of trajectory do not lie in one 
plane. Since a jet with a two-dimensional trajectory is a special case of the jet with the 
more general three-dimensional trajectory, it is desirable to have a theory that can esti- 
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mate the more general situation. At present, the only theoretical method which appears 
capable of predicting jet flows with three-dimensional trajectories is that of Hirst 
{ref, 21), although no results of this kind were presented in his report. 

In view of the preceding comments, the purpose of the present theoretical investi- 
gation is to develop an integral method which accounts for natural fluid mechanisms such 
as turbulence, entrainment, buoyancy, and heat transfer in the conservation equations 
governing the jet flow. In particular, it is desirable to have a theory that (1) utilizes the 
momentum conservation equation along the jet path, in addition to the one perpendicular 
to the path, in order to avoid any assumption regarding the growth of the jet cross- 
sectional area as was done in references 10 and 15 to 17; (2) obtains a third momentum 
conservation equation for the jet flow which, when solved simultaneously with the other 
two momentum equations, allows three-dimensional jet paths to be calculated, thus extend- 
ing the theories reported in references 18 and 22; (3) provides estimates of the tempera- 
ture of the jet fluid by examining several heat-transfer mechanisms that can account for 
the heat loss from the jet flow; and (4) is easily adapted to account for free -stream flows 
with either a nonuniform velocity field or a nonuniform temperature field. 

THEORETICAL DEVELOPMENT 

This section is concerned with the development of a theory which approximates the 
fluid mechanical process that occurs when a turbulent jet of circular cross section is 
injected into a semi-infinite free-stream flow. The mathematical model allows the jet to 
penetrate into the cross flow and to bend over and spread under the influence of natural 
fluidic forces, and the jet velocity vector to approach the free-stream velocity vector at 
some point downstream from the jet exit. This capability is obtained by considering a 
section of jet fluid as a control volume similar to the approach used by Reilly (ref. 16). 

The control volume is illustrated in figure 2 which depicts the trajectory of a jet injecting 

into a free-stream flow having a velocity taken 

to be spatially but not time dependent. The origin of 
the Cartesian (x,y,z) coordinate system is at the injec- 
tion point, whereas the natural (s,n,t) coordinate sys- 
tem moves and rotates as it follows the path of the jet 
axis which is traced out by the jet velocity vector. The 
s-axis is located along the trajectory, the n-axis is 
oriented perpendicular to the trajectory in the direc- 
tion of the radius of curvature of the trajectory, and 
the t-axis is perpendicular to both the s- and n-axes 
and completes the right-hand axis system. A complete 
discussion of the natural coordinate system is given in 



Figxire 2.- Cartesian and natural 
coordinate systems. 
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appendix A. The equations expressing conservation of mass and momentum in the direc- 
tion of the natural axe&are derived in the following sections. 

Continuity 

The following integral expression equates the net influx of mass into the control 
volume to the rate of increase in mass in the control volume: 

where d6 and do represent the elemental surface area and volume, respectively. It 
is assumed that the flow process is steady (in the mean), fully turbulent, and incompress- 
ible. The fact that incompressibility is assumed does not imply that the flow process is 
one of constant density. Equation (1) thus becomes 

Carrying out the operations suggested by this equation leads to the mass flows through the 
surfaces of the control volume CD ^)' 

It is noted that d6 = dA for surfaces © and ©, which represent a cross section 
of the jet flow, that is, areas that are perpendicular to the trajectory. Since mass flow 
is a continuous single- value function of position along the trajectory, a Taylor expansion 
can be performed to obtain the mass flow through surface © as a function of th^ass 
flow through surface ©. The difference between mass flows through surfaces ® and 
(© represents the mass flow through the sloping face ©, which defines the amount of 
free- stream fluid mass (per unit of jet length) that is entrained, or drawn, into the 
control volume. Entrainment can be written as 


E 



p V dA 
I ^ 


(3) 


The two most common profiles used in the literature to describe the velocity varia- 
tion at a given cross section in the jet flow are the Gaussian and ”top-hat." The Gaussian 
representation is particularly useful for providing estimates of flow properties on the jet 
center line and was used in the development of Hirst (ref. 21). In using the Gaussian pro- 
files however, the assumption is made that the jet flow is circular, that is, axisymmetric, 
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which is valid for free-jet flows and for coflowing jet flows but is questionable for a jet 
injecting into a cross flow where the jet cross section is not circular (ref. 6). In addi- 
tion, the theoretical results obtained by using a Gaussian velocity profile are only appli- 
cable in the region where the jet flow has become fully developed. 

Top-hat profiles represent the average jet flow properties and have been used in 
a number of theoretical studies (refs. 10, 15, 16, 18, and 20). By using the average flow 
properties in the conservation equations, it is not necessary to place restrictions on the 
symmetry of, or to assume similarity in, the jet flow. This means that the governing 
equations can be used to describe the jet flow in the region where the jet is fully devel- 
oped or in the region where the jet is only partially developed and a potential core still 
exists. Accordingly, the present study uses the averaged jet flow properties defined as 
follows: 



which state that at a given location along the trajectory the jet local velocity, density, and 
temperature values are integrated over the jet cross-sectional area. These definitions 

permit the conservation of mass in the jet (eq. (3)) to be expressed in differential form 
as 


^ _ d(pAV) 

" ~~dT~ (5) 

Continuity equations of this form are obtained in references 12, 16, and 18. 

In recent years it has been recognized that inclusion of the entrainment process in 
a jet injection analysis is important not only because of its influence on jet momentum, 
and hence trajectory, but also because this process allows for the mixing of jet and free- 
stream scalar properties such as temperature. The fact that entrainment occurs when 
there is relative motion between two flow fields has been used by different researchers 
to justify relating entrainment to the appropriate velocities normal and parallel to the 
jet axis. This type of representation of the entrainment function E results in a variety 
of empirical -constants” which must be adjusted in order to obtain suitable agreement 
with experimental data. Examples of several attempts to approximate E in this fashion 
are seen in references 8, 19, and 20. The model of entrainment presented in reference 18 
uses three constants, one of which is obtained by satisfying the Ricou- Spalding measure- 
ments (ref. 13) for an isothermal free jet. Experiment trajectory data, rather than mass 
flux information, were used as the criterion for indirectly adjusting the other two con- 
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stants. Experimental mass flow data are presented by Kamotani and Greber in refer- 
ence 8 to support the constants used in their model of the entrainment function. 

Because of the complex helical streamline pattern (usually interpreted as a pair of 
counterrotating vortices) evident in the lee side of the jet, it is believed that the entrain 
ment function cannot be split into two totally independent parts as described in the previ- 
ous paragraph. One attempt to account for the free-stream fluid entrained into the jet as 
a result of the helical circulation pattern has been reported by Flatten and Keffer (ref. 23) 
Their entrainment model was extended by Hirst (ref. 21) to account for local buoyancy in 

the jet flow. 

In consideration of the complicated nature of analytically predicting the entrainment 
function and the desire to keep empirical constants to a minimum, the present study 
defines this jet property by using the experimental data for an air jet obtained in refer- 
ence 6. The entrainment function is given by 



P„E*(V - VJ 


( 6 ) 


where C is the effective jet circumference. Equation (6) represents the entrainment 
parameter in functional form, where the entrainment coefficient E* was obtained from 
measurements of mass flux along the jet axis (ref. 6). In order to use this formulation 
for E, V must always be equal to or greater than The entrainment coefficient is 

presented in figure 3 for three injection velocities and, as noted, is a function of both 
injection velocity and distance along the trajectory. The effective velocity ratio Vgff, 
defined as the square root of the ratio of injection to free-stream dynamic pressure, has 


been suggested in references 3 and 8 as the 
proper parameter with which to compare dif- 
ferent injection situations, and hence is uti- 
lized throughout the present paper. An 
empirical expression was obtained here to 
represent E * in the present mathematical 
model and is shown compared with the exper- 
imental data. It is assumed that this func- 
tion can be used at larger values of Veff and 
s ydj than those shown in the figure. 

It should be noted that the empirical 
expression for E* implies that as s — 0, 

E * (and hence E) - 0. This is an unset- 
tling possibility because, near s = 0, the jet 
flow in this type of injection process might 


Experiment (Keffer and Baines (ref. 6)^ 



Figure 5-- Entrainment coefficient as a 
function of position along axis of 
isothermal air jet. = 9^^* 
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be expected to resemble that of a free jet with similar entrainment rates. Keffer and 
Baines (ref. 6) observed that as s - 0, E* was of the same order as that found by 
Morton (ref. 12) for simple jet flows. This fact is further substantiated if the work of 
Ricou and Spalding (ref. 13) who measured entrainment for axisymmetric turbulent free- 
jets is considered. The formula they suggested as best representing their data is 


= 0.32 |- 

(paV)i di 


(7) 


This equation can be differentiated with respect to s to get the entrainment in the free 
jet. Putting the results in the same form used by Morton gives 


E = O.OBCiPjVi ,8) 

Thus, the constant 0.08 compares favorably with Morton's value of 0.116. 

It is assumed in the present study that E* = 0.08 until the empirical expression 
presented in figure 3 predicts a greater value, which is then used. This is reasonable in 
view of the fact that an increase in Vgff decreases E* so that as E* 

should approach the free- jet value. 


Conservation of Momentum 

The following integral expression equates the rate of increase of momentum in the 
control volume to the sum of forces acting on the control volume plus the net influx of 
momentum into the control volume: 


I iTj” = HI ‘b + jj ? • - II V(P V . (9) 

wher^the surface stress tensor t is taken to contain both shear and pressure terms 
and f |3 = Body force per unit volume. Since a time independent flow process has been 
assumed, the term on the left-hand side of equation (9) is zero; hence, 

III 'b Od + II “ • - II . e,^)de - 0 (10) 


This equation will be used to obtain expressions for momentum in the direction of the 
natural coordinates. 


U 



Assuming the body forces are independent of the integral over do gives 




( 11 ) 


where, as a first approximation, the volume is assumed to be equal to A As, which neg- 
lects the rate of change of A with s. The length of the jet control volume is repre- 
sented by As. Since the components of f^ and later are desired in the s-, n-, 

and t-directions, a vector identity can be used to obtain 

% " (^b ■ ®s)®s (^b ■ ®n)®n (^b ‘ ®t)®t 

and 

V« = (V«, • eg)es + (V^ • + (v^ • e^)e^ (1^ 


where 

There are two body forces considered to be acting on the jet flow. The first force 
is due to the buoyant condition which results from a difference in density between the jet 
and free-stream fluids. This force acts in the y-direction and is given by 


The second force results from centrifugal effects associated with the jet having mass and 
following a curved path. This force acts in the direction of the radius of curvature of 
the trajectory and is expressed as 


- ^ _ -pV2 

f - f e ss — — 


(15) 


The total body force per unit volume on the control volume f^^ is the vector sum of the 
buoyancy and centrifugal forces and is substituted into equation (12) to obtain 

?b = + [«g(by ■ f„) + ‘c]bn + 'gfy ' ?t)«l 

The pressure portion of the shear-stress tensor term in equation (10) can be writ- 
ten in component form as 


15 



(17) 


P/s '>« P,^s <•« +'0 P/n «« 

® @ ® 

® ® 

where p^ represents the local pressure on the respective surfaces indicated by the 
integration. 

n-momentum .- The n-niomentum equation is obtained by taking the n-components 
of the various vector quantities in the momentum integral equation (10). The term repre- 
senting the body forces, for example, is extracted from equation (16) as 

^b,n = ^g(®y ■ ®n) + (18) 

and the pressure force is obtained from equation (17) to be 

^P.n ° iT P^ (IP) 

® 

This pressure force is combined with the shear stress integrated over surface to 
obtain the total drag force on the jet flow due to blockage of the free- stream flow. 

It is noted that for many injection situations the centrifugal force is in the opposite direc- 
tion from the drag force, the drag force being in the positive n-direction. 

Because of the complexity of the interaction between the jet and free -stream flows, 
the force resulting from the blockage effect is postulated to be the drag on an equivalent 
solid cylindrical jet shape inclined at an angle to the free- stream flow (after Abramovich 
(ref. 10)). This force can be expressed as 

°n = CD,n^„^„S^ef,n (20) 

where the dynamic pressure of the free- stream flow perpendicular to the jet axis is given 
by 


q 


°o,n 


1 

2 


P V 


2 

po,n 


Employing the definition • e„) results in 


(21) 
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(22) 


q 


00 ,n 



The reference area for calculating the drag force is defined by 


®ref,n 


= h AS 


(23) 


where h is the local width of the jet measured in the t-direction. Incorporating the 

expressions for ^ and into equation (20) allows the drag force to be finally 

expressed as 


D„ =C 


^n 2 


D,n^ 


o(^x • ^n) 


h As 


(24) 


The third term in equation (10) accounts for the net influx of momentum into the 
control volume, which for the n-direction, is due solely to the flux across the slanted 
surface of the cylinder. It is represented by the rate at which mass enters the sides 
of the control volume (eq. (3)) multiplied by the free-stream velocity component in the 
n-direction. Thus, 



% d6) = 


,Uv 

j oo,n as[J J 


P,V,dA’ 


As = V E As 
«>,n 


(25) 


Integrating equation (18) over the control volume and combining the results with equa- 
tions (24) and (25) yield the following conservation of momentum equation in the n-direction: 


pAv2 

R 



- P)(5y • "n)A + Cp „q_^h(e^ • e„)2 , EV„(e^ ■ r„) 


(26) 


The expressions for the dot products in terms of the Cartesian coordinates as the depend- 
ent variables and the distance along the trajectory as the independent variable are given 
in appendix A. Substituting the appropriate definitions the n-momentum equation becomes 


PAV2 

R 


d^y 


= gA(p^-P)R^+Ci,„q 




+ EV R 


d^x 

ds^ 


(27) 


where the radius of curvature R of the trajectory is defined in appendix A. Previous 
studies of the 2-D injection problem (usually vertical injection) proceed to put this equa- 
tion into a form with the slope of the trajectory, or the angular orientation, as the depend- 
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ent variable. Since the vertical and lateral injection problems are treated herein as spe- 
cial cases of the more general 3-D injection situation, it is desirable not to specialize the 
equations in that way. The reader will see in appendix B how this equation is nondimen- 
sionalized and put into the form used in obtaining a numerical solution on the computer. 
The direction cosines of the jet velocity vector (i.e., were chosen as the dependent 
variables for several reasons, one of which is that the order of the governing equations 
is reduced by one. In addition, the algebra is thereby kept to a manageable level. 

It is appropriate to mention that at this stage of the development, the approach used 
in a number of previous studies (refs. 10 and 15 to 17) has been to assume the area growth 
of the jet along the trajectory, the rate of growth being based on data measurements where 
s/dj 5 10, This was done in lieu of solving the s-momentum equation. The area growth 
can be obtained by assigning a certain shape for the jet cross section (e.g., circular or 
elliptical) and by allowing the jet width to grow at a specified rate. Schetz and Billig 
(ref. 15) used the expression for mass flow in the jet, that is, the continuity equation, to 
eliminate velocity in the n-momentum equation, the resulting expression then being inte- 
grated to obtain a solution for the jet trajectory. Typical results obtained by this pro- 
cedure are presented in figure 4 for a jet with an elliptical cross-sectional shape and are 
compared with experimental data acquired from the photographs in reference 17. 

Comparison of the assumed cross-sectional areas with the values obtained from 
experiments shows that the two are in reasonable agreement in the proximity of the jet 
exit but that the values diverge as the jet proceeds downstream, the measured areas indi- 
cating a much more rapid rate of jet growth than the assumed values. This trend is 
reflected in the trajectory information where good agreement between the predicted and 
experimental trajectories is noted in the initial region after jet injection, but poorer 
agreement occurs farther downstream. It should be mentioned that the investigators 
who made the area-growth assumption were predominantly interested in the jet trajec- 
tory in the proximity of the injection point. Another effect of assuming area growth is 
seen in the erroneous trends for the theoretical jet velocity, illustrated in figure 4 by the 
velocity deficit curve. As noted, the jet velocity begins to increase at some point along 
the trajectory. The reader is aware, of course, that as long as the jet injection velocity 
IS greater than the free-stream velocity, the jet velocity will decrease continuously along 
the trajectory and eventually approach the free-stream velocity value far downstream. 

From these remarks it is obvious that an alternative approach should be considered 
such that the jet cross-sectional area is permitted to be an unknown in the governing equa- 
tions. In order to do this it is necessary to have another equation to solve along with the 
continuity and n-momentum equations. The equation expressing conservation of momen- 
tum along the trajectory satisfies this need. By using this additional momentum equation, 
a more natural description of the jet flow properties is obtained as illustrated in figure 4. 
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- Experimental and theoretical Jet flow prop- 
for. a water injection process4 = 18 . 2 ; 

= 1.0; - 90°- 


Although the area assumption is to be discarded, it is still necessary to provide 
information concerning the width of the jet in order to calculate the drag terms in the 
governing equations. One approach is to follow Abramovich (ref. 10) and assume the 
growth of the jet width along the trajectory by using an empirical expression based on 
limited experimental data. Another approach is to specify a shape for the jet cross sec- 
tion and to use this with the computed area to calculate the jet width. This latter approach 
is more appealing because it is easier to justify its use on the basis of available experi- 
mental data. Keffer and Baines (ref. 6), for example, have shown that a jet initially hav- 
ing a circular cross section transforms to a '^kidney*^ shape as the jet penetrates into 
the cross flow. This shape remains approximately the same with increase in s as 
illustrated in figure 5. 

Prior attempts by researchers at approximating the jet cross section have been 
limited to the elliptical shape, where the circle is a degenerate case. Hirst (ref. 21), 
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Figure 5 .- Cross-sectional pressure contouirs of jet injected, 
normally into cross flow, V^/V^ 2.2 (after Abramovich 

(ref. 10)). Solid and dashed lines represent lines of 
constant total and static pressure, respectively; shaded 
areas indicate potential core region. 

for example, approximated the jet cross section as a circle, which was a useful assump- 
tion in his development because then he assumed the jet flow to be axisymmetric. It 
would appear that the ellipse would be a more suitable approximation, particularly near 
the injection point where the jet flow is deformed by the large pressure and shear stress 
fields. An ellipse with a major-to-minor axis ratio of 5:1 was employed in reference 10 
and later in references 15 to 17, whereas reference 18 assumed a 4:1 ellipse. The approx 
imation used in reference 18 had the added advantage of accounting for the change from 
the circular shape at the jet exit to the 4:1 ellipse at a specified point along the trajectory. 

For the present study, the jet cross-section shape is assumed to be elliptical with 
a ratio of major-to-minor axes of 5:1, the major axis being the jet width. It is also 
assumed that at the injection point (s = 0) the elliptical area is equal to that of a circle 
with diameter dj^. The change of jet width h along the trajectory is accounted for by 
the expression 



The change in h with distance along the trajectory was obtained by using calculated val- 
ues of A resulting from the solution of the governing conservation equations. This vari- 
ation IS shown in figure 6 compared with experimental data obtained from Kamotani and 
Greber (ref. 8). The data points from reference 8 were obtained from contours of con- 
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Figure 6.- Variation of jet width with distance along trajectory for a 
range of injection velocities, = 9^^, “^iAa. " 

stant velocity normal to a jet cross section, where the velocity excess had decreased to 
10 percent of the maximum excess. The theory predicts the trend of increase in h/d^ 
with increase in s /di, although the predicted values are higher than the measured values 
at large s/dj. The small effect of Veff on h/d^ near the injection point is reflected 
by the theory. The empirical expression used by Abramovich (ref. 10) in his theory is 
shown for comparison. 

The value of Cj^ associated with the 5:1 elliptical shape is taken to be 1.6 in 
keeping with the equivalent "solid" body argument and is assumed to be independent of 
the Reynolds number of the flow over the ellipse Wooler, Burghart, and Gallagher 

(ref. 18) used a value of 1.8 in their analysis, whereas Abramovich (ref. 10) used 3.0, a 
value which was pointed out by Schetz and Billig (ref. 15) as being totally unrealistic. An 
indication of how sensitive the theory is to the choice of cross-section shape and blockage 
coefficient can be seen in figure 7 where experimental and theoretical trajectory results 
are compared for = 4.0. The effect of changing ellipse axis ratio is presented in 

figure 7(a) where Cj) n = theoretical calculations, whereas the effect 

of varying Cj^ n specific ellipse axis ratio (5:1) is shown in figure 7(b). As noted, 

increasing either ellipse axis ratio or blockage coefficient results in progressively lower 
theoretical trajectories. The trends discussed here for Vgjf = 4.0 are typical of other 
injection velocities. The fact that the theoretical results obtained for Cj) ^ =1-® 

(fig. 7(b)) are in better agreement with the experimental trajectories than the results 
obtained for Cj) „ = 1.6 is misleading. It will be demonstrated later (fig. 22) that the 
theory estimates mass flows in the jet that are too low, so that if this disparity was 
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(a) - 1,6. (b) ellipse. 


Figure 7*- Effect of jet cross-sectional shape and blockage coefficient on 
theoretical trajectories for =4.0, = 90 *^ j = 0°. 

corrected the results obtained with Cp = 1,6 would show improved agreement with 
the experiment. ’ 

s-momentum .- The s-momentum equation is obtained by taking the s-components 
of the various vector quantities in equation (10). The resulting expression represents a 
balance between the rate of change of jet momentum and the forces on the jet due to 
changes in mean jet pressure, to buoyancy caused by density differences between the 
jet and free-stream fluids and to entrainment of ambient fluid into the jet. 

In the s-direction the shear and pressure forces are not combined as they were in 
the n-direction. The appropriate pressure force acquired from equation (17) is 

j'i'p,''' (29) 

® (D CD 

where the local pressures are integrated over the control surface to get 
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( 30 ) 


A Taylor expansion can be performed to obtain the pressure and area at surface ® in 
terms of the appropriate variables at surface (T), whereas taken to be the 

average of pq and p^ and Aq to be the difference between Aq and A0 

Substituting these expressions into equation (30), neglecting terms having higher orders 
of ds, gives the pressure force in the s-direction as 


The force contribution of the surface stress tensor in the s-direction is zero because the 
local jet velocity at the control volume surface is equal to the free-stream velocity com- 
ponent tangent to the jet flow, thus negating a shearing stress. 

The contribution of the net influx of momentum into the control volume to the force 
balance in the s-direction is obtained by performing the operations suggested by the third 
term in equation (10) and is found to be 




ds) = ^(pAV2)as - EV 



(32) 


The definition of averaged jet properties was used to derive this equation, and the momen- 
tum flux entering the sloping surface of the control volume was represented by the rate 
that mass flows across the surface multipled by the free-stream velocity component in 
the s-direction. Equations (31) and (32) are combined with the body force term from 
equation (16) and the dot product expressions from appendix A to obtain the following 
s-momentum equation: 


9(pAv2) 




— + EV — 
as ds 


(33) 


In order to evaluate the static pressure gradient along the trajectory (ap/as), the 
assumption is made that the free-stream static-pressure field around the jet perimeter 
imposes itself on the jet flow. This is the usual type of assumption made concerning 
other free turbulent processes, such as jet injection parallel to a mainstream (coflowing 
flow) or jet injection into a reservoir (free jet flow). For the present case where the jet 
structure is considered, as an elliptical cylinder inclined at an angle to the free-stream 
flow, there are large variations in the free-stream pressure field around the jet due to 
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Figure 8.- Static pressure variation around 
the perimeter of a circular cylinder. 


the blockage effect that the jet has on the 
free-stream flow. Some idea of the static 
pressure variation around the perimeter of 
a jet cross section, idealized as a circular 
cylinder, can be obtained by observing the 
experimental pressures in figure 8. (These 
experimental pressures were obtained from 
ref. 24.) An assumed pressure distribution 
to be used in the theory is also presented. 
As noted, the assumed pressures on the 
front of the cylinder (o i 0 s are in func- 
tional form and were obtained from poten- 
tial flow theory, whereas the pressures on 
the back of the cylinder ^2 ^ ^ 
assumed to be equal to the free-stream 
pressure. Several researchers (e.g., 


Ramsey (ref. 7) and Kamotani and Greber 
(ref. 8)) have approximated the pressure 
field around the jet by examining the potential flow over various cylindrical shapes. 
Although the pressure field resulting from the turbulent jet injection process is very 
complicated and does not lend itself to be categorized in this simple a fashion, the pres- 
sure variation in figure 8 is adequate for use in the present mathematical model. The 


local surface pressure (Cpq^ ^ + p J is used in the expression 


P = 



(34) 


to obtain the average static pressure acting on the cylinder. Performing the integrations 
in equation (34) with the pressure distribution shown in figure 8 gives 


where it is recalled that q^^^ is the dynamic pressure resulting from the free-stream 
velocity component normal to the trajectory. (See eq. (22).) This equation implies that 
the average static pressure on the jet cross section is less than the free-stream static 
pressure but approaches p^ as q^ ^ approaches zero. This occurs when 
approaches zero and/or when the jet Becomes parallel to the free-stream flow (i.e., 
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a = 0). If the average pressure is assumed to impose itself on the jet flow (i.e., the 
pressure in the jet flow becomes p), p can be differentiated with respect to s to get 


ds 


= -<1. 




ds 


Vds^ 


+ R 


2 d^x d^x 
ds^ ds^ 


(36) 


It is noted that if the jet trajectory is in either the vertical (X-Y) or lateral (X-Z) planes, 
the expression for the pressure gradient (eq. (36)) simplifies to 


ds 


da 

sin 0 ! cos a — 
ds 


(37) 


which is the form used by Campbell and Schetz (ref. 22). Incorporating the pressure gra- 
dient term (eq. (36)) into equation (33) yields the following final form of the s-momentum 
equation: 


3(pAv2) 


dy 


as 




R 


dR/Vx^ 


ds 


\ds' 


ds^ ds3 


EV^ p- 
ds 


(38) 


t-momentum.- The t-momentum equation is obtained by taking the t-components of 
the various vector quantities in equation (10). The resulting expression represents a bal- 
ance of forces on the control volume due to buoyancy, to blockage of the free-stream 
flow, and to entrainment of ambient fluid into the jet. Similar to the derivation of the 
n-momentum equation, the pressure force is combined with the shear stress integrated 
over surface to obtain the total force acting on the jet in the t-direction. 
Accordingly, 

•^t ~ ^D,t‘l=«,t®ref,t 


where 



n the jet cross section were circular, then Cq ^ would equal and Sj.gf would 

equal However, since the shape is elliptic and not circular, is not equal 

to Cjj and t ®ref,n’ elliptical shape has an axis ratio 

of 5:1, Sref ,t = I ^s, so that 


25 



( 41 ) 


where Cq j. is taken to be 1.0. 

The net flux of momentum in the t-direction entering through the sloping surface of 
the control volume is represented by the rate at which mass flows across the surface 
multiplied by the free-stream velocity component in the t-direction; this results in 

't ■ il • fi)As (42) 

The body force term from equation (16) is combined with equations (41) and (42) to obtain 
the t- momentum equation: 


gA(p^ - 


• ®t) * ^D,t 1 


e^)^ + EV 


' ®t) = 0 


(43) 


After substituting for the dot products (appendix A), the torsion associated with the 
trajectory is arranged so that it is in the numerator of the terms. The reason for this 
arrangement is that is expected to have a small value for the present study Xt = 0 
for a two-dimensional trajectory The resulting expression is ^ 


S ^ ^ f (P„ - P)S 


d^\ 


gA, 




+ Cn,q +Cr. q llR2/'d3sf^P 

+ i l4 hR 3^^ ^+-Cr,.qhSS^^+ic a !l3eA<1x 


+ T EV ^ ^ + T EV R & + EV la dx _ 


o °° ds j„2 ' ’0 ' «>-■ 3 


ds' 


ds"- 


R ds 


(44) 


For the case where the jet follows a two-dimensional path, the t-momentum equation is 
an identity (see appendix B); hence, its use is not necessary in the procedure for obtain- 
ing a solution of the jet trajectory and flow properties. 


26 



Heat Energy 


Until this point only the mass and momentum aspects of the jet injection process 
have been discussed; however, since the present investigation is concerned with heated 
discharges, it is necessary to also consider appropriate methods of describing the ther- 
mal characteristics of the flow. In particular, it is advantageous to determine the change 
in mean jet temperature resulting from the penetration of the jet into the cross flow. This 
determination can be accomplished by monitoring the heat loss from the control volume, 
the heat loss resulting from several heat-transfer mechanisms. 

The first type of heat-transfer mechanism pertains to the reduction in energy con- 
tent per unit volume pCpT of the jet fluid due to the entrainment of free-stream fluid at 
a different energy level (PCpT) ■ Applying this concept to the control volume results in 
the expression °° 

(mCpT)^ = {“‘=pT)i ^ 

where (mCpTj represents the energy level in the control volume that would exist if 
there were no entrainment and (mCpT)^ represents the equilibrium energy level result- 
ing from the complete mixing of the jet and entrained fluids. The various specific heats 
in equation (45) are assumed to have the same value. 

Forced convection, the second type of heat-transfer mechanism being considered, 
results when the free stream flows around the heated jet fluid and extracts heat energy 
from the jet in the process. This heat transfer is analogous to the forced convection in 
separated flow over a heated cylinder, where the cylinder is cooled by the fluid flowing 
normal to the cylinder axis. To be consistent with our previous approach, the convective 
heat transfer is estimated by considering the jet structure as a cylinder inclined at an 
angle to the free-stream flow. 

Eckert and Drake (ref. 25) give several examples of film heat-transfer coefficients 
occurring in this type of flow situation and suggest the following expression for estimating 
an average Nusselt number: 

%u,d = 

The value of Prandtl number for an air injection process is taken to be 1.0, whereas for 
a water injection process the functional dependence of Np^, on water temperature was 
obtained from tables in reference 25. The Reynolds number is defined with the "effective 
diameter of the jet as the reference length and the free-stream velocity component per- 
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pendicular to the jet axis as the reference velocity so that Nusselt number will be sensi- 
tive to the changes in local flow conditions as the jet penetrates into the cross flow. The 
Reynolds number is thus 


N 


Re,d 


. PooVoo.nd 




V ,2 
oo ds 


For vertical or lateral injection this equation reduces to 


"Re,d “ V 


a = sin a 


(47) 


(48) 


The values of occurring in the present study suggest the selection of :// = 0,45 

and I = 0.50 for use in equation (46). 

The definition of Nusselt number is Hd/k where k represents the thermal con- 
ductivity of the jet fluid. This definition is used to obtain the average film heat-transfer 
•coefficient H which yields the following rate of heat loss from the jet fluid: 


Q = HS^(T^ 



(49) 


where denotes the cylindrical area of the jet control volume. This, in turn, results 
in a temperature change in the jet flow due to this convective heat loss. 

An example of the temperature results obtained when these two heat-transfer mech- 
anisms are incorporated into the solution of the governing conservation equations is shown 
in figure 9, where the theoretical calculations were made with the same injection condi- 



Pisure 9*- Variation of Jet temperature with 
distance along trajectory for heated 
injection process. > T„. 


tions as the water injection process reported 
in reference 22. The trend of temperature 
decrease along the trajectory measured in 
reference 22 is adequately estimated by the 
theory, the predicted average temperature 
values falling below the measured maximum 
temperatures. These results are substan- 
tiated by the temperature data for air injec- 
tion processes measured by Ramsey (ref. 7) 
and Kamotani and Greber (ref. 8). The theo- 
retical temperatures obtained by considering 
only the effects of entrainment are presented 
to demonstrate the relative magnitudes of the 
two types of heat-transfer mechanisms. Con- 
vection is seen to be the dominant mechanism 
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for determining temperature loss in the early stages of the jet injection process, whereas 
the effects of entrainment become dominant as the jet proceeds downstream. 


Solution Procedure 

This section summarizes the highly nonlinear governing differential equations and 
briefly describes the iterative method that is employed to obtain a numerical solution at 
specific locations along the jet trajectory. Appendix B shows how the equations are non 
dimensionalized and put into the forms used in the numerical technique. The governing 
equations are repeated here for convenience as follows: 

Continuity (eq. (6)): 

E=4p^E*(V-V„) 


n-momentum (eq. (27)): 




d^x 

ds2 


5 -momentum (eq. (38)): 


a(pAv2) xdy . A 


r, dR /d^x\ . o2 d^x 
ds2dp3 


+ EV^ ^ 
ds 


t-momentum (eq. (44)): 
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^ d^ + T EV R + EV ^ ^ = 0 
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Heat transfer due to 


entrainment (eq. (45)): 

(mCpT)^ = (mCpT)^ + ”a('pT)^ 

convection (eq. (49)): 

Q = HSp(T„ - T) 


It was found that the s- , n- , and t-momentum equations could be simplified some- 
v/hat by using direction cosines u and w as the dependent variables rather than x, 
y, and z. This procedure means that at each point j+1 on the trajectory a solution to 


this initial value problem involves determining values for u- i, w- i o A anri 

]+!’ ]+!’ j+1’ “““ 


The basic solution procedure is to solve the s-momentum equation for the jet 
momentum in the control volume, where the coefficients in that equation are estimated 
by using the flow property values obtained from the solution at the previous location on 
the trajectory. The s-momentum equation is used in conjunction with the continuity equa- 
tion to provide an update on and the heat loss from the control volume is calculated 

to provide new estimates for and The current flow property values are then 

used in the coefficients of the n-momentum equation to obtain a solution for (du/ds)- 
where a central difference scheme provides Uj^g- deration between the s- and 



n-momentum equations provides the information 
used to solve the t-momentum equations for 
(d2w/ds2).^j, from which wj^g obtained by 
using a central-difference scheme. The most 
recently calculated values of the direction cosines 
and flow properties are used to iterate back through 
the governing equations. (See flow chart of itera- 
tion process in fig, 10.) For the case where the jet 
path is two dimensional the information from the 
t-momentum equation is redundant; therefore, it 
is only necessary to solve the s- and n-momentum 
equations , along with the continuity equation, to 
obtain the desired solution. It was found that con- 


’Number Of fteratfons 

Figure 10.- Flow chart illustrating 
iteration oF governing equations 
at a specific location on jet 
trajectory. 


vergence to a satisfactory solution occurred in only 
a few iterations for a two-dimensional trajectory 
case and in less than 10 for a three-dimensional 
trajectory case. 
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Incremental values of x, y, and z are obtained from the final value of trajectory 
slope and from the assigned value for As. These increments are added to the coordinates 
of the previous location on the trajectory to obtain new x,y,z trajectory coordinates. 

This procedure is repeated at each incremental "step" along the trajectory to provide a 
solution for the trajectory and cross-sectional area of the jet, as well as the jet flow 
properties of mass, velocity, momentum, and temperature. The theoretical results pre 
sented herein were obtained with a constant incremental step size of O.Old^. 

EXAMINATION OF THEORETICAL RESULTS 

The purpose of this section is to discuss some of the limitations of the theory devel- 
oped in the previous section and to demonstrate its versatility for handling a variety of 
injection situations. In order to establish the authenticity of the present theoretical 
method for estimating jet flow properties, its predictions are compared with estimates 
from other analytical models as well as with experimental data acquired from a number 
of studies. The last portion of the section presents a theoretical example of a jet with a 
three-dimensional trajectory. 


Two-Dimensional Trajectory 


Experimental trajectory data obtained 
from different investigations of air jets are 
presented in figures 11 and 12 and show the 
two-dimensional paths of the turbulent jet 
for a range of injection velocities and ori- 
entations. These data were obtained from 
hot-wire measurements and, thus, represent 
the path that is traced by the maximum 
velocity in the jet flow. The measurements 
show that an increase in injection velocity 
or angle Qfj results in farther penetration 
of the jet into the free-stream flow. Theo- 
retical trajectories were calculated with 
the same injection conditions as the exper- 
imental data and are in good agreement with 
the measured trajectories throughout the 
range of injection velocities and orientations. 

Experimental trajectories are pre- 
sented in figures 13, 14, and 15 for a water 
injection process where the jet is injected 
perpendicular (ai = 90°) to the direction of 



Figure 11.- Experimental and fbeoretical 
trajectories of an air Jet having a 
range of injection velocities. 

= 90°; Pi = 0°. 
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Figijre 15.- Experimental and theoretical Fi^e l4.- Experimental and theoretical 

trajectories for lateral injection pro- trajectories for oblique injection 

cess. ATi = -1.1 K; = 90°; = 90°. process, = 90°; 3i = 50°. 


flow in a water channel. The trajectories were measured from photographs of the jet 
injecting laterally (0^ = 90°), obliquely = 50°), and vertically = 0°) into the main 
stream, and, thus, they represent the path of the approximate center line of the jet flow. 
The information for the oblique injection process (fig. 14) is presented as the projection 
of the trajectories onto the vertical (X-Y) and horizontal (X-Z) planes. The effect of 
increasing injection velocity on the trajectories of these water jets is the same as that 
observed for the air jets; that is, the jet penetrates farther into the main stream with 
increase in injection velocity. 

Theoretical trajectories calculated with the same injection conditions as the experi- 
mental data adequately represent the trend of increasing jet penetration with injection 
velocity and are in good agreement with the measured trajectories for the lateral and 
oblique injection processes (figs, 13 and 14, respectively). However, the present theory 
(represented in fig. 15 as a solid line) does not predict the amount of penetration experi- 
enced by the jet injected vertically. The reason for this discrepancy is that the experi- 
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mental results of figures 13 and 14 were obtained for completely turbulent jet flows, 
whereas the data in figure 15 were acquired for jet flows that are laminar at the injec- 
tion point. It is recalled, of course, that the theory was developed by assuming fully 
turbulent flow. In order to estimate the trajectories for this mixed-flow situation, the 
present theory was adjusted to account for the initial laminar portion of the jet flow. 

This adjustment was accomplished by assuming that the jet begins its turbulent growth 
at a point y /\ specified in the photographs in reference 26. Since the location and 
extent of the transition region in the flow are functions of injection conditions ^e.g., 

^Re di}» well as free-stream conditions (e.g., it is expected that the values 

of ^y A change accordingly. The appropriate values of y^A modify 

the theory are shown in figure 15 and the resulting calculations represented by the 

dashed lines. 

Theoretical trajectories calculated with the present theory are compared with theo- 
retical and experimental results of other researchers in figures 16 to 19. The analytical 
methods of Abramovich (ref. 10), Schetz and Billig (ref. 15), and Reilly (ref. 16) provide 
jet trajectories which are comparable to those of the present theory for the injection con- 
ditions presented (figs. 16 and 17). It is recalled that the theories of references 10, 15, 
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Figi^e 16.- Comparison of theoretical tra- Figure IT.- Comparison of theoretical trajec- 
^eetories estimated by present theory tories estimated by present theory with 

and by Abramovich (ref, 10). = 90 °; those of Schetz and Billig (ref. I 5 ) and 

^i ^ • Reilly (ref. 16) . = 90°; = 0°. 

and 16 assumed the growth of the jet cross-sectional area along the trajectory by using an 
empirical expression (fig, 6) based on experimental data for s/d^ ^ 10 . As a conse- 
quence, the trajectories predicted by these theories agree quite well with experimental 
data in the vicinity of the injection point. Care must be exercised in using these theories 

to estimate jet trajectories and flow properties at large s/dj^ values. 

One of the best theoretical methods prior to the present is that of Hirst (ref. 21), 
who attempts to account for the complex flow processes that take place as the flow evolves 
from a momentum jet near the injection point to a buoyant plume at large distances down- 
stream, His results are compared with the present theory in figures 18 and 19 for a 
range of injection velocities and angles. As noted, the present theory is in better agree- 
ment with the bulk of experimental data for all of the Injection conditions. Since Hirst 
assumed a Gaussian type of velocity distribution in the jet, his theory is applicable only 
in the region where the jet flow has become fully developed. This explains why his theo- 
retical trajectories do not originate at the injection point. The experimental data obtained 
by Gordier (ref. 2) are shown in figure 18 because Hirst compared his theory with these 
data in reference 21. Gordier's data, however, indicate greater penetration by the jet 
than is seen for the other data. Ramsey (ref. 7) suggested that this discrepancy was 

probably due to injection into a cross flow with a very thick boundary layer. This trend 
is shown later. 
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Figure l8.- Comparison of theoretical trajectories 
estimated by present theory with those of Hirst 
for normal injection, ct^ = 90° j Pi = 0°. 
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Jet Flow Properties 

Examples of sonte of the theoretical flow properties obtained in the process of solv- 
ing the governing conservation equations are presented in figures 20 to 25. Experimental 
jet areas are shown In figure 20 normalized by the jet area at the injection point and 
plotted as functions of s/d^. These meas- 
urements indicate that the cross-sectional 
area of the jet continually increases as the 
jet proceeds along the trajectory. The data 
points acquired from Kamotani and Gerber's 
work (ref. 8) were obtained by measuring 
the area encompassed by a contour of jet 
velocity where the velocity excess had 
decreased to 10 percent of the maximum 
excess with respect to the free-stream 
velocity component tangent to the jet flow. 

It is noted that Keffer and Baines (ref. 6) 
also defined the edge of the jet flow in 
this fashion. Experimental areas for 
Veff nQUj = 30 were obtained from the pho- 
tographic information in reference 22 by 
assuming the cross-sectional shape to be a 
5:1 ellipse and by measuring the minor axis. 

The data in figure 20 also indicate that an 
increase in injection velocity results in 
larger rates of area growth with s/dj. 

Theoretical areas are shown for comparison and predict the same trends as the 
experimental areas; however, the theory underestimates the magnitude of the area growth 
experienced by the jet for the range of injection velocities shown. It should be mentioned 
that these theoretical estimates of jet area are very sensitive to the amount of entrain- 
ment, small increases in entrainment resulting in large increases in the cross-sectional 
areas. It will be shown shortly from mass-flow arguments that the entrainment model 
used in the theory could be improved, which would result in better agreement between the 
theory and the experiment in figure 20. 

As the area occupied by the jet fluid grows with increase in distance along the tra- 
jectory, the jet velocity correspondingly decays; this is illustrated in figure 21 where the 
jet velocity is nondimens ionalized by the injection velocity and where the V^/V^ values 
for the various injection conditions are depicted by the arrows. The experimental data 
that are shown are measurements of the maximum jet velocity for Vgff ^om = ^ 


Theory Experiment 

Veff, non 

_ 0 4 8 



Figure 20.- Variation of jet cross-sectional 
area with distance along trajectory for a 
range of injection velocities. 

(Values from ref. 8 are measurements of 
the area bounded by a contour of jet 
velocity) . 
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and for a free jet. These data indicate that 
the trend for velocity decay is similar for 
all injection velocities. 

For all the cases, there is a short dis- 
tance (potential core) where the maximum jet 
velocity remains equal to the Injection veloc- 
ity. Increase in s/dj beyond this point 
results in continued decreases in jet velocity 
which eventually approaches the free-stream 
velocity value (note arrows in 

c 10 20 30 40 fig. 21). Increasing injection velocity 

increases the potential core length, a maxi- 

Fi^re 21.- Effect of injection velocity on mum value being obtained for the free iet 

jet velocity decay with increase in dis- / . 

tance along trajectory, - 90 ^* (ref. 8), and decreases the value of ^oo/^t 

that the jet velocity must approach, 

Voo/Vi = 0 for the free jet. The combination of these experimental trends helps in 
comprehending the theoretical velocity decay curves presented in the figure for a range 
of injection velocities. Comparing the theoretical velocities with the experimental data 
for Vgff jjQjjj = 4 and 8 shows that the theory predicts y/Vj^ values that are greater 
than the experimental values and that approach at a slower rate. The fact that 

the mathematical model estimates average jet velocities explains why the theory predicts 
a decrease in velocity immediately after injection compared with the existence of a poten- 
tial core demonstrated by the maximum velocity measurements. 

The effects of increased injection velocity on velocity decay, which were discussed 
in connection with figure 21, are put into perspective in figure 22 by plotting the velocity 
as a velocity deficit, where the difference between the jet and free-stream velocity is 
divided by the difference at the injection point. Presenting the results in this fashion 
causes the velocity deficit to approach zero as V ^ Y^. As noted, an increase in 
^eff,nom results in a corresponding increase in experimental velocity deficit at a given 
distance on the trajectory. This trend also applies essentially to the theoretical velocity 
deficit variations. Keffer and Baines (ref. 6) observed that velocity deficit showed a uni- 
versality [i.e., independent of Vgff] when plotted against the distance from the virtual 
source of the jet flow. 

The variation of theoretical mass flow in the jet with distance along the trajectory 
is shown in figure 23 for a range of injection velocities. As would be expected from our 
consideration of mass conservation in figure 3, the mass flow in the jet increases with 
increased distance along the trajectory, the higher mass flows occurring for the lower 
injection velocities. This increase in mass flow with increase in s^d^ is demonstrated 



40 




Figure 22.- Variation of jet velocity 
deficit with distance along tra- 
jectory for range of injection 
velocities, = $ 0 ° . 


Theory 


Experiment 
eff.nom 



Schetz and Billfg (ref. 15) 
/ 

t i 

40 


Figure 23 .- Variation of mss flow in jet 
with distance along trajectory for range 
of injection velocities, = 90 ®. 


by experimental data (ref. 8) for Vgff^nom = although the measured values indicate 
higher mass flows in the jet than are predicted by the theory. A further indication that 
the theoretical mass flows should be greater than they are Is provided by the measure- 
ments of mass flow in a free jet. (See ref. 13.) The fact that these mass flows are less 
than the experimental values for the jet with Vgff = 4 supports the premise that the 
entrainment rate for a jet in a cross flow should be greater than that for a free jet. At 
certain values of s/dj the disturbing situation exists that the theory predicts mass flows 
that are less than the free-jet values. Improvements in the entrainment model used in 
the current analytical effort will provide improved estimates of the jet mass flow and 
consequently will yield more realistic cross-sectional area and velocity decay results 
than were observed in the last few figures. In the study by Schetz and Billig (ref. 15), 
the mass flow was assumed to remain constant at the Initial value. Their assumption 
is represented by the horizontal dashed line in the figure and becomes more realistic 
nearer the injection point. 

The variation of theoretical jet momentum with distance along the trajectory is 
shown in figure 24 for a range of injection velocities. The jet momentum, normalized 
by the value at the injection point, is seen to decrease immediately after injection and 
reach a point on the trajectory where it attains a minimum value, after which it increases 
for the remainder of the trajectory. This trend occurs for a jet injected normal to a free- 
stream flow because the component of momentum in the y-direction initially decreases 
along the trajectory faster than the component in the x-direction increases. For this 
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situation = 90° and the y- and s-momenta are identical at the injection point, so that 
the natural decrease of y-momentum with increase in s^dj results in a decrease in 
s-momentum during the initial phase of the injection process. As the jet axis becomes 
parallel to the free -stream direction, the y-momentum of the jet approaches zero and the 
X- and s-momenta become synonymous. 

The continual decrease of the y-momentum along the trajectory of a jet injected 
normal to a cross flow is illustrated in figure 25, where experimental data of Kamotani 
and Greber {ref. 8) are presented for several injection velocities. As you would expect, 
the y-momentum is largest ^at a given s/dj location) for the jet with the highest injec- 
tion velocity. These trends of y-momentum with increase in s^di and V^ff are also 
reflected by the theoretical results of the present study. Abramovich (ref. 10) assumed 
in his analytical development that the component of jet momentum perpendicular to the 
free-stream direction (y-momentum) remains constant along the trajectory. The fallacy 
of this assumption is particularly obvious at large distances. 



Figure 24 .- Variation of theoretical jet 
momentum with distance along trajec- 
tory for range of injection veloci- 
ties. = 0 ^. 



Figure 25.- Effect of injection velocity 
on the change of jet y-momentum with 
distance along trajectory, a. = 90^: 
Pi = 0°. 


Free-Stream Nonuniformities 

One advantage of the present theory is its flexibility for Investigating parameters 
which affect the trajectory and flow properties of the injected jet. Not the least important 
of these parameters are the free-stream velocity and temperature fields into which the 
jet is injected. Up until now the free-stream velocity and temperature have been assumed 
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constant, but some of the effects resulting from relaxing these assumptions are demon- 
strated in figures 26 to 29. A water injection process is used for this demonstration. 

For the purpose of this illustration, the jet is assumed to inject vertically into 
a free-stream flow which has a boundary-layer type of velocity distribution in the 
y-direction (fig. 26(a)). A Karman-Pohlhausen velocity function is described from the 
injection surface to the boundary- layer edge, which is taken to be 8dj^. At larger val- 
ues of y/dj the velocity is assumed to be constant having the same value as used 

for the uniform free-stream velocity case shown in figure 26(b). Figure 27 shows that 
injection into the nonuniform free-stream velocity field results in farther penetration by 
the jet into the cross flow than injection into the free stream with the uniform velocity 
field. Coincident with this, the jet velocity decay is essentially unaffected, whereas the 
jet cross-sectional area and momentum are less at any given distance along the trajec- 
tory. It was noted that the effect of free-stream velocity nonuniformity on jet trajectory 
diminished with increase in injection velocity. 



(b) Uniform free-stream velocity. Figure 27.- Effect of free-stream velocity 

nonuniformity on theoretical trajec-^ 

Figure 26.- Uniformity of free-stream veloc- tory. = 0^. 

ity field, = 90°; = 0°. 


The linear temperature gradient shown in figure 28 is used to demonstrate the 
effect of injecting a heated jet vertically into a free-stream flow having constant velocity 
and nonuniform temperature fields. The free-stream temperature at the injection sur- 
face (295.4 K) is equivalent to the value used for the uniform free-stream temperature 
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case. Injecting a jet with initial temperatures of (295.4 K) and (306.5 K) into a free- 
stream flow having a uniform temperature results in the trajectory and temperature 
curves shown in figure 29. The combination of the jet flow having a very large Froude 
number and experiencing a rapid heat loss results in no change in the trajectory. This 
result agrees with the experimental observations made by Kamotani and Greber (ref. 8). 
As noted in the figure, the temperature for the heated jet decreases along the trajectory 
until it reaches the free-stream value Injection of the heated jet into the free- 

stream flow with a temperature gradient results in a trajectory similar to that obtained 
by injecting the heated jet into the uniform temperature field. There is a definite differ- 
ence, however, in the temperature curves resulting for these two injection conditions, 
where the nonuniform T^ situation results in higher jet temperatures because of the 
larger values of T^ that the jet flow "sees” as it penetrates into the cross flow. Heat 
is initially lost from this jet flow until a point is reached on the trajectory where a heat 
gain is experienced. 




x/d| 


(b) Uniform free-stream temperature. 




^Figure 28.- Uniformity of free-stream tern- figure 29*- Effect of free-stream temperature 
perature field, = 90 ^; ^ 0 ^ . nonuniformity on theoretical jet trajectory 

and temperature. Veff =5; ^i = 90 ° i 
Pi = 0°. 
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Three-Dimensional Trajectory 

Up until now all the theoretical trajectories that have been presented are, by defini- 
tion, two dimensional; that is, they lie in a single plane. Figures 30 and 31 are presented 
to illustrate a three-dimensional path that results for a water injection process where a 
jet is injected normal (a^ = 90°) to the mainstream and rotated 45° away from the vertical 
An injection situation similar to this was observed in figure 14, where the theoretical tra- 
jectories are actually two dimensional because the Froude number of the jet flow 
large; or in other words, the buoyancy force is small compared to the momentum forces. 
In order to obtain a three-dimensional trajectory, particularly near the Injection point, 
the Froude number of the jet flow must be small. 


For the purpose of this demonstration, the heat transfer from the jet is ignored so 
that the jet density is assumed constant along the trajectory. Solutions for the governing 

conservation equations were obtained for Npr,i Pi = Poo) %r,i " 

the projections of the resulting trajectories on the X-Y and X-Z planes are presented in 
figure 30. For i - projections on the two planes are equivalent which occurs 

only when /Sj = 45° "'and the trajectory is two dimensional. Allowing to be less than 
p leads to the second set of projections, 
which shows that decreasing the Froude num- 
ber increases the penetration of the jet into 
the cross flow. The fact that the effect is 
more pronounced on the X-Y projection than 
on the X-Z projection is an indication that 
the triad of unit vectors associated with the 
natural coordinate system is tracing out a 
three-dimensional path. In effect these unit 
vectors are "twisting" out of the plane where 
the two-dimensional trajectory obtained for 

— oo is located. 

^ Fr,i 

Further evidence that these numerical 
results are consistent can be obtained by 
examining the dependent variables of the 
governing equations (u and w) and their 
gradients. This information is presented 
in figure 31 where the injection conditions 
are identical to those used in figure 30, and 
the parameters in question are plotted as 



functions of distance along the trajectory. 
The solid curves represent the case where 


Figure 50.- Effect of Froude number on jet 
trajectory for =4. = 90° j 

3i = 45°.. 
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Figiire 31*- Effect of Fronde number on variation of jet trajectory parameters 
with distance along trajectory for - 4. 3^ = 45°, 

buoyancy force is zero typical of the trends that result for injections 

with 0° S Q!j i 90° and 0° £ j3^ g 90°. As noted, an increase in s/dj^ results in a con- 
tinual increase in u and a decrease in w, u and w approaching 1,0 and 0.0, respec- 
tively, for large s/dj. These limits coincide with the jet velocity vector which becomes 
parallel to the free-stream velocity vector. In conjunction with these trends for u and 
w, du/ds and dw/ds are, respectively, positive and negative valued. 

The dashed curves indicate what happens to the trajectory parameters when a siza- 
ble buoyancy force is considered. The fact that the values of u and w which result 
when = 10 are, respectively, smaller and larger than the corresponding values 

when Npj. j is indicative of the increase in penetration experienced by the jet when 
the buoyancy force is added. An interesting aspect of this situation is seen in the varia- 
tion of w with s/di, where w reaches a minimum value at some point on the trajec- 
tory and then begins to increase. This is the same as saying that achieves a maxi- 
mum value and then begins to decrease. The trend of w with s/dj is reflected in the 
dw/ds curve, where dw/ds changes from its usual negative sign to a positive value. 
This result is important because it signifies that buoyancy force is becoming dominant 
over the other forces acting on the jet flow. It is intuitively obvious that for the situation 
where buoyancy is the driving force, the jet path will tend toward the vertical so that w 
will be increasing with increase in s/dj and dw/ds will be greater than zero. This 
can be shown explicitly by letting = u = du/ds =0 in the n-momentum equation pre- 
sented in appendix B {eq, (B9)) and by solving the resulting expression for dw/ds. 
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It is not possible to substantiate the trends discussed herein with results from other 
studies, since no experiinental data exists for a jet with a three-dimensional trajectory. 

In the theoretical development of reference 21, it was stated that the governing equations 
could apply to a jet following a three-dimensional path; however, no calculations of a 
three-dimensional trajectory were presented to support that claim. 

CONCLUDING REMARKS 

An investigation has been conducted to develop a theoretical method which accounts 
for the pertinent fluid mechanic and heat-transfer aspects of heated turbulent jets dis- 
charging into a cross flow. The injection process was assumed to be incompressible, 
where the jet flow injects into a cross flow of semi-infinite extent having a spatially 
dependent velocity field. 

A theoretical model of the jet was developed by using an integral method, which 
accounts for natural fluid mechanisms such as turbulence, entrainment, buoyancy, and 
heat transfer. Solving the governing conservation equations yielded predictions of jet 
trajectory and area growth that agreed well with experimental results for air and water 
injection processes having a wide range of injection conditions. Because the present 
investigation obtained the jet cross-sectional area in the process of solving the governing 
equations (unlike previous studies), the present theory allowed a prediction for various 
jet flow properties, such as velocity and momentum, to be obtained and provided better 
estimates for the jet trajectory. 

Realistic estimates of temperature in the jet fluid were obtained by accounting for 
heat losses in the jet flow due to forced convection and to entrainment of free-stream 
fluid into the jet. Forced convection was seen to be the dominant heat-transfer mecha- 
nism during the early stages of the jet injection process, whereas the effects of entrain- 
ment became dominant as the jet penetrated further into the free-stream flow. 

The versatility of the theory was demonstrated by observing the effects of a jet 
injected into free-stream flows with either a nonuniform velocity field or a nonuntform 
temperature field. Theoretical results were also shown to illustrate a truly three- 
dimensional jet trajectory which was calculated by considering the Injection Froude 
number to be small. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 25, 1973. 
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APPENDIX A 


SPACE CURVE INFORMATION 


This appendix provides additional Information about the natural coordinate system 
used in the analytical study and presents the expressions necessary to transform the 
momentum equations in the text to the form used in the numerical solution. 


Space Curves 



(a) Natural coordinate system. 



Figure 32.- Illustration of the 
natural coordinate system 
and the direction cosines 
of the unit vector 


The theoretical development utilizes a natural 
coordinate system that is associated with the path of the 
jet flow which consists of a moving triad of unit vectors 
eg, e'j^, and ej.. This coordinate system is illustrated 
in figure 32(a) where r is the position vector from the 
origin to a point on the curve and is given by 

r=xe^ + yey + ze^ (Al) 

The derivative of a position vector is shown by Hildebrand 
(ref. 27) to be a unit vector tangent to the curve and point- 
ing in the direction of increasing arc length as follows: 


®s 


= ^ = 
ds ds X 



y 


+ ^e 
ds z 


(A2) 


It is noted that dr/dt is the velocity vector associated 
with a point moving with speed ds/df along the curve. 
Thus, 


V=^^ = Veg 

dt ds ® 


(A3) 


Since e’g is a unit vector, it follows that 



(A4) 


The derivative of with respect to s has a direction perpendicular to the 
curve and is written as 
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ds 'ds2 ^ 


C[2y ^ 

ds 


+ 

ds2 " 


(A5) 


where the length of this vector is the curvature of the curve. Defining the radius of cur 
vature R as the reciprocal of the curvature and e„ as the unit vector in a direction 
normal to eg gives 


^ ^ L r (A6) 

dsR^ 


so that 



(A7) 


A moving, rotating triad of mutually orthogonal unit vectors is completely described 
by the addition of the third unit vector e^ which by definition is 



X e 


n 


(A8) 


It is noted that for a plane curve, and lie in the plane of the curve, whereas e^. 
is a constant unit vector perpendicular to that plane. The trajectory resulting from the 
vertical injection process is an example of this situation. Differentiating equation (A8) 
leads to 


del _ 
ds"*' "Vn 


(A9) 


where the scalar is the torsion of the curve, the negative sign implying that is 

positive when the vector triad rotates in a right-handed sense about eg as it progresses 
along the curve . 

In order to find de^j/ds, write x Og and differentiate to get 



^o®t - 


R s 


(A 10) 


Equations (A6), (A9), and (AlO) are known as the Frenet-Serret formulas. Taking the dot 
product of ej. with den/ds leads to an expression for torsion which can be written in 
determinant form as 
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dx/ds 

dy/ds 

dz/ds 

d^x/ds^ 

d2yy^s2 

d^z/ds^ 

d^x/ds^ 

d^y/(ls3 

d^z/ds^ 


(All) 


The dot products that are required in the governing momentum equations can now 
be found. For example, e^ • e„ and e^ • are needed in the n-momentum expression 
(eq. (26)). Using equations (A 5) and (A 6) yields 


e„ = R 


d^r 

ds^ 


(A12) 


so that 


and 




(A13) 


(A14) 


A similar procedure is used to get the dot products needed for the s-momentum 
expression; 


and 



dx 

ds 



(A15) 


(A16) 


The process for obtaining the dot products for the t-momentum expression (eq. (43)) 
is more involved because of the use of e^.. The equation for ej. is found from equa- 
tion (AlO) to be 


(A17) 
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Taking the dot product of and Cy, respectively, with this expression and using the 
distributive law for dot products leads to 


^ ^ ^ ^ ^ + -i— ^ (A 18 ) 

t Tq ds ^g2 Tq ^g3 RTq ds 


- - _ 1 J_^ (A19) 

■ ®t - To ds ds2 ^o ds3 ds 


Direction Cosines 

It was mentioned in the section "Theoretical Development" that the procedure for 

solving the governing equations (6), (27), (38), and (44) simultaneously was simplified by 

using the direction cosines of the unit vector eg as the dependent variables. The angles 

e makes with the X-, Y-, and Z-axes are shown in figure 32(b) and are defined as 
S 


COS a 


dx 

ds 


cos, = | 


cos 



(A20) 


Equation (A4) represents the auxiliary equation for these direction cosines and is used to 
express the third direction cosine (dz/ds) in terms of the other two. Define 




so that equation (A4) becomes 

— = (l - u^ - w2) = 77^/2 

ds 


The derivatives of this expression are found to be 


(A21) 


(A22) 


d2z 

ds2 


w^ 

ds ds _ 


V 


1/2 


A 

„l/2 

V 


(A23) 


and 
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-u 



.ds/ ds2 



(A24) 


Expressions for R, dR/ds, Tq 

At this point it is desirable to obtain the expressions for R, dR/ds, and in 
terms of the direction cosines u and w. Equations (A21) and (A23) are substituted into 
equation (A7) to get the follov/ing expression for R: 


1 

R 



+ 




(A25) 


The relation for dR/ds is obtained by differentiating R (eq. (A 7)) and substituting the 
necessary auxiliary equations from equations (A21) to (A24) as follows: 


-L ^ ^ d^u . ^ d^w 
r 3 ds ds ^jg2 ds ds2 


, 1/2 {, 1/2 


dg2 \ds/ ds2 




(A26) 


The torsion of the curve is obtained after some manipulation by expanding the determinant 
in equation (All) and by substituting the auxiliary equations as follows: 









+ - 


V 


A 

1/2 


VL + 7} 


1/2 /_±_ ^ 1/2 A 

\ O ^ 1 /O \ O 


ks^ 77^2 


ds 


ds 


(A27) 


Thus, equations (A25) to (A27) represent R, dR/ds, and Tq as functions of u and w 
and their derivatives for the general case of a three-dimensional curve. It is interesting 
to note the forms these equations take when the trajectory is two dimensional, such as for 
a vertical injection process (y = 90°). For this situation, the auxiliary equation (A22) 
becomes 


u^ + w 


2 = 1 


(t? = 0) (A28) 
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so that 

A^u^+w^ = 0 <A29) 

ds ds 

Using these relations in equation (A25) leads to 

(A30) 

R w ds 

Likewise, 

1 dR _ du d^u dw d V (A31) 

■ r 3 ds " ds ^jg2 ds ds2 

and Tq = 0. Torsion would be expected to equal zero from equation (A9) since the tra- 
jectory is a plane curve which implies that e^ is a constant vector (parallel to the 
Z-axis). 
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NONDIMENSIONAL CONSERVATION EQUATIONS 

The purpose of this appendix is to nondimensionalize the governing conservation 
equations and to make several observations concerning their numerical solution. In 
addition, special forms of the governing equations are examined and the t-momentum 
equation is shown to be an identity for the vertical injection situation. 


s- Momentum 

The direction cosine expressions from appendix A are substituted into equation (38) 

(pAV^)t 

and the resulting expression is divided by — ^ to get the nondimensionalized 
s- momentum equation as follows: 


1_ d(pAV2) _ 1 M - P)” ^ 1 /A VU f = dg/duf ^ p2 du A 

y2\ ds 2 \Aij ^ q^ 2(Aijl^qjy|_ ds W 


(pAv2\ 


+ j(^ u) (Bl) 


where the barred symbols indicate division by dj. The expression 


E = J 


(pAV)i 

di 


(B2) 


was also used in the process of obtaining equation (Bl). 

As was mentioned in the section ’’Solution Procedure,” equation (Bl) is solved at 
point j+1 for the jet momentum which is estimated by using a backward finite difference 
to approximate the rate of change of pAV^ with s. This gives 


d(pAv2) 


ds 


_ (pAV^)j,l - (PAV^)] 


'j+1 


As 


(B3) 


which can be substituted into equation (Bl), The momentum is cal- 

culated by using values for the various parameters in equation (Bl) evaluated at j as a 
first estimate, and then using j+1 values after the first iteration. The jet momentum 
is then used in conjunction with (pAV)j^j to provide values for Vj^j and, hence, Aj^j. 

At this point in the numerical solution the heat loss from the jet control volume is 
accounted for. It should be recalled (see eq. (45)) that the mass of fluid in the control 
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volume must be specified before and after entrainment takes place in order to calculate 
the effect of entrainm^t on jet temperature. The total mass in the control volume after 
entrainment is given by 

"^3+1 = + 

where the entrained mass is found to be 

_ (AS)^E (B5) 

“^e " V 

The velocity of entrainment Vg is taken to be an average velocity in the control volume, 
equation (6) is used in conjunction with the E* function presented in figure 3 to specify 
E. 


n-Momentum 

The expressions for the direction cosines (appendix A) are substituted into equa- 
tion (27) and the resulting equation is divided by pAV^. The ratio of mass flow in the 
jet to the initial jet mass flow is written as 

K - (B6) 

(pAV)i 

and is used to get the nondomensional n-momentum equation as follows: 


L = g,r^+g,rY^'| +g,r^ 

A ds 


R 


\ds 


ds 


(B7) 


where 




/P.V A f ~ P) 


^ 2h r f^(R.\(A.\ 

^ ' ,k2 “'VlAPliW 


Gs 
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In order to put equation (B7) into the form used in the numerical solution, the 
expression for R must be considered (eq. (A25)). Examining the vertical injection case 
results in a cubic equation in du/ds after equations (A25), (A28), (A29), and (A30) are 
substituted into equation (B7). It is noted that two roots of the cubic equation are zero if 
the buoyancy term is neglected. This type of trivial solution can be avoided by dividing 
du/ds out of equation (A25) to get 


1 

R 





1 - u2 - w2 


1/2 



(B8) 


which is substituted into equation (B7) to get the following transformed n-momentum 
equation: 


du _ Qj dw ^3 ^4 

ds . 2 du .3 A 2 
Aj Aj Ai 


(B9) 


This equation is used to calculate (du/ds).^j, with w.^j, (dw/du)j^j, Gj, G2, 

and Gg being known. Appropriate values of u are then obtained by utilizing 

(1) A backward finite -difference scheme at the first step away from the injection 
point (i.e., s = As) 

“i-i ' “i * (BIO) 

(2) A central finite-difference scheme to get Uj^g points on the remainder of 
the trajectory (i.e., s > As) 

= “i * (Bll) 

For the vertical or lateral injection situations, the auxiliary relation (eq. (A22)) can 
be used to obtain a value of wj^2 corresponding to the value of Uj^2 equation (Bll), 
while its derivative 


u^ + w ^ = 0 
ds ds 


(B12) 
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can be used to obtain a value of (dw/ds)j^j and, hence, (dw/du)^^j^. These parameters 
are used to iterate through the governing equations to obtain a solution for the two- 
dimensional trajectory. For the more general injection cases, the t-momentum equa- 
tion is required to provide information on w and its derivatives. 

t- Momentum 

The direction cosine expressions from appendix A are substituted into equation (44) 
and the resulting relation is divided by d^/q^ to obtain the nondimensional t-momentum 
equation as follows: 


FT ^+FoT ^+F„T +F4f^f + F,-fe) +^6 + ^7^^ 
o ds 2^0 ^-2 3 o ^ 4\ds / 6 7 ds ds^ 

+ Fg ^ + Fg II + FiqTq II + F^^To ^ + F 12 T 0 = 0 
ds 


(B13) 


where Tq = r^dj^ and 


F - ^fo - o\^ 
" d^q^roo ^)ds 


■'2 = - '>) 


r, = i^(p -p) 
3 diq^R\ / 


IT -hr roo\(dVi 

^4 " 5 ^D,tl q^ l(^ds 


V li r 


'^o.\u2 




TP 2h p du 


Tp 2tl 

Fg - 


_ _ 2h r /^°o\u ^ 

9 " 5 iJr ds 


_ 7T J^oo\d\i 

^10 "2 W/df 


/V \ 

-p - 2L Ti p 

^11 ■ 2 Ui/ 
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It has been seen in equations (A26) and (A27) that dR/ds and are functions of 
d^w/ds^ and d^u/ds^ plus lower order terms. The approach used here is to define 
d^n/ds^ by using a finite difference approximation after a solution has been obtained 
from equation (B9). The problem then becomes one of getting dR/ds and Tq in terms 
of d^w/ds^. Thus, from equation (A26), 




and from equation (A 2 7), 


= Bq + B, 


/d2 


w 


\ds" 


(B14) 


(B15) 


where the B coefficients are functions of u, w, du/ds, dw/ds, and d^u/ds^. When 
these relations are substituted into the t-momentum equation (eq. (B13)), a quadratic 
expression in d^w/ds^ results as follows; 



+ B 


d^w 


+ = 0 


(B16) 


The coefficients are calculated from information obtained from the n-momentum equation, 
and a solution is acquired at point j+1 by approximating d^v//ds^ with a central finite 
difference. This results in a quadratic in which is solved to provide an update on 

Wj+2’ {dw/ds)j^l, and (dw/du)-_^j. 

It is desirable to demonstrate that the t-momentum equation is an identity for the 
vertical injection situation where the trajectory is two dimensional and is confined to the 
X-Y plane. For this case, e^ is parallel to the Z-axis which implies that the dot prod- 
ucts of this vector with e^^ and Cy in equation (43) are zero; hence, the t-momentum 
equation is an identity. This fact can be shown in a more rigorous fashion by considering 
the expressions for ^e^ • and • ej-^ from equations (A18) and (A19). Substituting for 
direction cosines into these two equations gives 


d^u dR du u ? fx 
2 ds ds R “ ^ 


(B17) 


and 
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n d^w . ^ ^ 4 . w ? ft (B18) 

ds2 dsds-^R- 

where the desire is to prove that these relations are equal to zero for the vertical injec- 
tion case. To do this it is necessary to use the auxiliary expressions that are presented 
in equation (A22), with ?j = 0, and in equation (B12). It is noted that equation (B17) can 
be written as 



(B19) 


The expression for R for a vertical injection case was shown in appendix A (equa- 
tion (A30)) and is substituted into equation (B19) to get 


du dw ? ft 
u ^ + w ^ = 0 

ds ds 


which is definitely zero from equation (B12). 

The procedure is similar for equation (B18) where 


(B20) 


d_ 

ds 




(B21) 


Equation (B12) is solved for dw/ds which, along with the relationship for R, is sub- 
stituted into equation (B21) to obtain 


_ ^ ^ = 0 (B22) 

ds ds 

Thus, the dot products e“jj • e^, and Cy • e!j. are zero for the vertical injection situation; 
hence, the t-momentum equation (eq. (43)) is an identity. This fact was used as a check 
of the computer results to see that the numerical output from the t-momentum equation 
was correct. 
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